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0^ , Abstract. Uncertainty relations for mixed quantum states (precisely, purity-boun- 

ded position- momentum relations, developed by Bastiaans and then by Man'ko and 
Dodonov) are studied in general multi-dimensional case. An expression for family of 
^ ' mixed states at the lower bound of uncertainty relation is obtained. It is shown, that 

in case of entropy-bounded uncertainty relations, lower-bound state is thermal, and a 
transition from one-dimensional problem to multi-dimensional one is trivial. Results 
of numerical calculation of the relation lower bound for different types of generalized 
purity are presented. Analytical expressions for general purity-bounded relations for 
highly mixed states are obtained. PACS number: 03.65.Ca. 



1. Introduction and review 



qh| Well-known position-momentum uncertainty relation for standard deviations of x and 

p operators, 

X ■ AxAp ^ h/2, (1) 

5r ! is valid for any state (described either by a wavefunction or by a density matrix [1,2]) and 
plays an quite an important role in quantum physics. In particular, uncertainty relation 
sets the precision limits of measurement process for non-commuting observables [3,4]. 
Another important example is that generalized coherent states (and squeezed states) 
could be defined as a set of states which minimize an uncertainty relation, see [5]. 
Uncertainty principle and properties of its minimum is also of special interest in theory of 
operators in Hilbert space, see further references in a recent work by Goh and Miccelli [6]. 

The inequality have been generalized to include extra dependence on degree of 
purity [7] of a quantum state 

^^ = Tr(p2) (2) 

(p is a density operator), the parameter ^ yU ^ 1 and equality // = 1 is achieved only 
for pure states. An asymptotic inequality for one-dimensional highly mixed states with 
/i ^ 1 has a form [8-14] 

AxAp>--. (3) 
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In addition to the trace Q of squared density operator, there are other measures of 
overall purity (see above cited papers for details, especially a recent comprehensive 
review on purity-bounded relations [14]). 

Another approach for treatment of uncertainty relation for mixed states are 
developed, by Wolf, Ponomarenko, Agarwal [15,16], and also by Vourdas and his co- 
authors [17,18]. In the cited works, the uncertainty relation is expressed in terms of 
correlations of respective observables. On the other hand, the inequality of the type Q 
relates uncertainties in conjugated variables and a measure of overall purity of state. 

A generalization of the uncertainty relation to multidimensional space (vector 
observables, which can appear e. g. for multimode states, or multi-particle situations) 
was investigated in early days of quantum mechanics [19] (see also a review in [13]) and 
is still drawing attention of researches [20-22]. In its most simple form, uncertainty 
relation for n-dimensional position and momentum operators X = {xi,X2, ■ ■ ■ ,Xn), 
P = (Pi,P25 • • • iVn) could be written as 

(AXAP)"^fjV, (4) 



with definitions 

^ n 1 " 

(Ax)2 = - rr Ax„ (AP)2 = - n a^,. 

n J- J- n J- J- 



n ^ ^ n 

1=1 i=l 



In fact, due to equality between different coordinates in minimum of uncertainty relation, 
the inequality (jH) has the same meaning as n-th degree of standard one- dimensional 
relation ((H), see also discussion in [23]. 

The problem of generalization of the inequality Q for multidimensional case was 
treated in papers by Karelin and Lazaruk [23,24]. In particular, in our first paper on 
this topic [23], it has been shown with the help of the Wigner function formalism, that 
there is a non-trivial dependence of the purity-bounded uncertainty relation limit on 
the number of dimensions. For highly mixed states with /i <^ 1, 

where the parameter C{n) characterizes distance from a minimum (/j,/2)" for pure states. 
In deriving (0), we assumed, that the Wigner function of the minimum-uncertainty state 
is nonnegative. 

In another paper [24], the structure of the density matrix near the lower bound of 
uncertainty relation was also found, using decomposition of the density matrix in terms 
of Fock states (which form an orthogonal basis with a minimal uncertainty) 

p = ^ am,m'|"^i) \m2) . . . |m„) (m'J (m'al . . . (m^|, (6) 

m,m' 

with the additional condition 



= Ictm.m'l^ = const. 
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where m = (mi, . . . ,m<j) is a "vector" index with integer nonnegative components. 

At the lower bound of uncertainty relation, density matrix in Fock representation 
is diagonal, am,m' = o,m,m5m,m'; coefficients am,m depend linearly on the 'norm' of 
vector index ||m|| = X]r=i"^«' number of Fock states in representation of p is 
finite. Coefficients am.m of this decomposition are degenerate, and their multiplicity is 
determined by norm ||m|| and dimensionality n: 

(n) ^ (||m|| +n-l)\ 
3\\m\\ niminifn- IV ^'^ 



The inequality obtained in [24] correctly describes the whole range of fi, including 
perfectly pure case fi = 1. In particular, in the interpolating form it becomes 

AXAP ^ ^ !^±1IM (8) 
2 n + 2 ^ ^ 

with the auxiliary real parameter L{fi) being a root of the transcendental equation 

(n + 2L)(n + l)!r(L) 
^ (n + 2)r(L + n + l) ' ^""^ 
where T{y) is Euler's gamma-function. 

It is also necessary to note that the inequality, mathematically practically the same 
as uncertainty relation, but with another physical meaning, is often used for classical 
wave fields, e. g. in optics [8-12,23]. Results of the present article, as well as of preceding 
papers [23,24] could be used, with appropriate change of notations, for classical partially 
coherent fields and sources (in 1-, 2- and 3-dimensional space [25]). 

The uncertainty relations (0), ((HI) could be further generalized in order to take 
into account the dependence of inequality minimum on eigenvalues of density operator. 
Preliminary report on this topic, with stress on partially coherent classical fields, was 
pubhshed in [26]. Obtaining such a relation, together with the study of its asymptotics, 
is the main aim of the present paper. 



2. Uncertainty relation for the diagonal representation of the density 
matrix 

Any density matrix p has a spectral decomposition [27] 

P = ^Pm\'ipm) {ipml, (10) 
m 

where pm are the eigenvalues, and lipm) are eigenvectors of the density operator, then, 
each of vectors \iprn) could be represented via outer products of one-dimensional Fock 
states I A;) 

\^m)= Yl 4:L..,.J^l)l^2)...|fc.), (11) 

where k^, i = 1,2, . . . ,n corresponds to ith one- dimensional subspace. 
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The right-hand side of the uncertainty relation (jD) is calculated using the method 
of papers [8, 13]. Core idea of calculation is the introduction of auxiliary observable 

= i [(AP) + ^(AX)2] , (12) 

which could be regarded as energy of some oscillator with unit frequency and mass "i?. 
The minimum of with respect to (for = AX/AP) is exactly left-hand side of 
uncertainty relation 

mm Em = AXAP. 

As far as Fock states are eigenstates of the harmonic oscillator with eigenvalue 2k + 1 
(in the one-dimensional case), then substitution of pO|) and pi|) into (fT^ leads to 

AXAP^^^^p^ J2 [2ik, + k2 + ...kr,)+n] 

m ki,k2,...,kn 

Now, due to isomorphism between the set of all positive integers and the set of 
combinations of n positive integers, it is possible to consider the coefficients ^^^L kn 
as elements of some unitary matrix {Akm}- Then, (in)equality (fT^ can be treated with 
lemma:}: from [10] to give 

h 1 

AXAP^-- Yl ['^ih + k2 + ...k 

n 

(14) 

k-i,k2,...,kn 

where eigenvalues of the density matrix are ordered in a non-increasing sequence. 

Dependence of the expression 2{ki + k2 + ■ ■ ■ + kn) + n on indices /ci, k2, ■ ■ ■ ,kn is 
degenerate: this expression takes the same values for several combinations of indices. 
Therefore it is possible to rewrite the (in)equality (fT^ as 

h 1 

AX AP ^ - - Y^i^k + n) Y Pm{k) (15) 

k m=0 

where the values Qm^ (degeneration multiplicity) are defined by the formula ((Tj), and the 
eigenvalues of the density matrix are collected in groups of gm^ terms. Expression (fT3j) is 
the main result of the paper, and it is the most general form of the uncertainty relation 
for mixed states (partially coherent fields) in a multidimensional space. This inequality 
relates a minimal uncertainty volume of a state to the spectrum of the density operator 
corresponding to this state. 



A 



(m) 

ki,k2,...,k„ 



.(13) 



3. Multidimensional purity-bounded relations 

Using the method from papers [10,11], it is possible to find a dependence of uncertainty 
relation limit on some characteristics of purity of quantum system. Usually, a family of 

t See also [14]; to be self-contained, the lemma is reproduced, together with necessary changes for 
multidimensional case, in [Appendix A| of this article. 
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"generalized purities" (Shatten p-norms or "generalized entropies" [22]), is used, which 
is defined as 

/.«=[Tr(p^/(^-^))]'^-\ (16) 

where r is an arbitrary (not necessary integer) real number with r > 1. Important 
special cases of /i*^^^ include /i*^^^ = /i ("usual" purity, see above), "superpurity" 



/iW = lim/i(^'), (17) 



when only the largest eigenvalue of the density matrix is taken into account, and also 
"entropy-based" purity degree 

/is = exp(-S'), S" = -Tr(plnp), (18) 

which is defined in terms of Shannon- von Neumann entropy S, and so leads to "entropy- 
bounded" uncertainty relation. As it is shown in [12], /is can be treated as a limiting 
case of definition (fT^ for r oo: = Hs- 

"Superpurity" and entropy-bounded uncertainty relations play a special role for 
one dimensional case: owing to continuous non-increasing dependence of fx^^^ on r for 
r > 1 [12], they are limiting cases of family of characteristics p6|) . It is also possible 
to show, that the non- increasing dependence of /i'-^-' on r remains valid in a general 
multidimensional case, see details in [Appendix B 



As it can be easily shown by Lagrange method, at the minimum of uncertainty 
relation, the definition (fT^ reduces to 

r-l 



/i 



(n)xr-/(T--l) 
k I 



(19) 



where 



ik — E Pm(fc) 



m=0 



and 



Y.^k = l. (20) 

k 

Then the (in) equality ()15|1 may be rewritten to 

h 1 
2 n 



AXAP^^^Y.^2k + n)ik (21) 

k 



and, in order to obtain the relation of the type (0) for given structure of the density 
matrix, (i. e. eigenvalues pm) it is necessary to find a minimum with respect to variables 

is. 

As the degeneration multiplicity ((Tj) has a rather complex form, the task of detailed 
study of uncertainty relation minimum has, in general, no analytical solution (the same 
as in one-dimensional case [14]). Besides interpolated and asymptotic inequalities, which 
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will be studied later in the article, it is possible to obtain analytical solution for the case 
of entropy-bounded relations. 

Using Lagrange method, it is easy to show, that a minimum of uncertainty product 
AXAP ()2H) for given entropy S is attained if the coefficients C,m are given by 

U = Ag^::^exp{-Pm), (22) 

where A is a normalization constant and parameter (3 depends on the entropy. Taking 
into account structure of the density operator at the minimum of general uncertainty 
relation [compare representation in the case of entropy-bounded relation the 

density matrix is taking the form 

pr = ^'i^...p^'. (23) 



n times 



Here 



p^P = J2e-^'\k){k\ (24) 

A:=0 

is a density matrix corresponding to minimum of one dimensional entropy-bounded 
relation. Parameter f3 can be found from solution of transcendental equation 

S = Pn '^P^^^ - n In (1 - exp(-/5)) , (25) 
1 - exp(/3) 

which is in accordance with appropriate equation for one-dimensional case [13]. 
Uncertainty relation then could be written as 

{AXAPr^PiS) (^y, (26) 



^2^ 

or, for highly-mixed states with S ^ 1 

(AX AP)" ^ exp(S) 0y (27) 

Obtained structure of eigenstate decomposition is factorized on solutions of one 
dimensional problem (see [11,13]), which leads to thermal state ()24p . In other words, 
the entropy-bounded uncertainty relation of position-momentum type has no additional 
effects for multidimensional cases. 

In order to study general low-purity case, it is possible to utilize the approach from 
Bastiaans' paper [10], which is based on generalization of the Holder inequality. The 
mathematical details are presented in [Appendix C| 

Let's define an "uncertainty function" C{jj,^'^\n) (see also (0)) as 

C(^('^),n) =/i(") (^i|2M + n-2 [/i('^) 5(M,n,r)]'^'}^ , (28) 

where 

BiM,ri,r)= Yl T^^^^?^(^-< (29) 
^ — ^ m — 1 ! ml 

0<m<M ^ ' 
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Figure 1. Uncertainty relation minimum for < 1, n = 
variants of the degree of purity: o's — r 1, x 's — r = 2, *'s 
(entropy-bounded relation). 



. 6 and different 
3, -l-'s — r — > cxD 



and additional real and positive minimization parameter M is introduced. 
1-d problem has known asymptotical solution [10] 

C{fi^'-\l) = 2[r/{r+l)Y, /i(^)<l. 

In the same way, highly mixed states could be treated analytically for arbitrary 
multidimensional far as limit of small fx^^^ requires M to be sufficiently 

large [14, 24] , it is possible to replace summation in formula ()29|1 by integration together 
with approximation of degeneration multiplicity by m"'~^/{n — 1)! with 

1 r^'^ 

B(M,n,r)^- -/ m^-^M - mV dm. (30) 

[n - 1)! Jo 

The last relation could be calculated analytically, 

/n+l \ -1 

B{M, n, r) ^ M"+"+^ JJ^^ + ^^^'^ 
Appendix D for details). Further minimization of relation (j^Hj) with respect to M, 



see 



after some tedious but quite elementary algebra gives an asymptotic variant of general 
purity bounded uncertainty relation 

(AXAPV > (I)" C(„,r) = + (32) 

which describes the whole range of n and r for /i <^ 1. Resulting dependencies of 
C{n,t) on r for n = 2,3 together with one-dimensional case are presented in figure 
and figure El It is seen, that the obtained expressions correctly describe whole range 
of parameter r, demonstrating decrease of uncertainty minimum with increase of r and 
leading to entropy-bounded relations at r ^ oo. 
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Figure 2. Uncertainty relation minimum for /i^''-' <C 1, n = 1, 2, 3 and different 
variants of the degree of purity (solid lines, r — 1,...,100). Dashed lines denote 
uncertainty minimum for entropy-bounded relation with r — > c» (for n = 1, 2, 3, 
rspectively). 



4. Concluding remarks 

To summarize, it is worth to note that two main forms of uncertainty principle (of 
position-momentum type) for mixed states in multidimensional space are obtained 
in the present paper. The first one (fT^ relates minimal uncertainty product with 
eigenspectrum of the density matrix and the second (fH^ is a general (asymptotic) 
purity-bounded uncertainty relation. In both cases, minimum of uncertainty product 
is obtained when eigenstates of the density operator are Fock states. In the case of 
purity-bounded relation, eigenspectrum of the density operator is defined by 

im 9t\M - mY'^ , O^m^M (33) 

Appendix C| In other words, spectral representation of the density matrix is a 



see 



finite sum of Fock states. Such state is definitely non-classical, see discussion in the 
Dodonov's paper [14]. Upon transition to r ^ cxd (entropy-bounded relations), the 
minimum-uncertainty state becomes thermal (j211), i- e. classical. More detailed study 
of minimum- uncertainty states structure will be subject of another publication. 

It is also necessary to note, that the results of the paper is applicable for analysis 
and characterization of entangled quantum states. Indeed, the spectral decomposition 
of density matrix is closely connected to the Schmidt decomposition of non-separable 
states, see, e. g. [28]. Approach to uncertainty principle for entangled states can be 
based on mathematically analogous case of uncertainty (reciprocity) relations for pulsed 
partially coherent classical beam [25]. 
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Appendix A. 

The lemma is reproduced here mainly for completeness of material, and in order to 
make the above analysis clearer. Initially it was presented in Appendix A of paper [10], 
it also could be found in [14]. According to [10], idea of this proof was initially proposed 
by M. L. J. Hautus. 

Let the sequence of numbers be defined by 

oo 

fem = ^ \amk\'^7k, (A.l) 
fc=0 

where 7o ^ 7i ^ • • • ^ 7a; ^ ■ ■ ■ and coefficients a^k satisfy the orthonormality condition 

oo 

^amkaik^^mh m,Z = 0, 1,... (A.2) 

fc=0 

One may consider the numbers bm for m = 0, 1, . . . , M as the diagonal entries of an 
(M + l)-square Hermitian matrix H = \\hij\\ with 

oo 

hij = ^ aika*f.-fk, i, i = 0, 1, . . . , M. (A.3) 

k=0 

Let the eigenvalues P oi H he ordered according to 

/3o ^ A ^ . . . ^ ^ . . . ^ /3m. (A.4) 

Prom Cauchy's inequalities for eigenvalues of a submatrix of a Hermitian matrix [29], 
we can conclude, that /3m ^ 7m (m = 0, 1, . . . , M) and hence 

MM MM 

^bm^^ hmm ^^Pm>^lm- (A.5) 

m=0 m=0 m=0 m=0 

Furthermore, with the numbers Am (or pm, in this article) satisfying the property 
Aq ^ Ai ^ . . . ^ A,„ ^ . . ., we can formulate the chain of relations 



M M M 



m— 1 



1=0 1=0 



Am&m — Ao&O + Am&m — Ao^O + K 

m=0 m=l m=l 

M m M-1 

= Xobo + X] Am ^ 6; - ^ Am+l ^ bi 
m=l 1=0 m=0 1=0 

M-1 m M M-1 m 

= ^ Ki'^bi + Xu'^bi - ^ Am+l k (A.6) 

m=0 1=0 1=0 m=0 1=0 

M M-1 m 

= Am ^ fez + ^ (Am - Am+l) ^ bi 

1=0 m=0 ;=o 

M M-1 m M 

^ Am 7i + X] (Am - Am+l) X] 7i = X] Am7r7 



Im- 

1=0 m=0 1=0 m=0 
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On choosing 7^ = 2n + 1 and taking the hmit M — 00, we arrive at the inequahty 
00 00 00 
^ A„ ^ \amn\^{2n + 1) ^ ^ A„(2m + 1), 

m=0 n=0 m=0 

which becomes an equahty if \amn\ = ^mn- 

In order to modify this proof to miltidimensional case, it is necessary to choose 
7m as 7ni = 2(mi + . . . + nis) + n (here m = (mi, . . . , m„) is a "vectorial" summation 
index), and then to take into account degeneracy of coefficients 7in. 

Appendix B. 

By analogy with Bastiaans' paper [12] for any r, g, with 1 < r < g, holds 



(r) 



r-1 



E (n)/)r/(r-l 
Urn m 

I 



Y.0 

m 

(g-l)/(r-l) 



m \ m 



r-1 



(r-g)/(r-l)- 



-1 r-1 



(B.l) 



Here 6^^ = ^m/gm\ and the Holder inequality for weighted sum [30] is used, see 
formula ()(y.3jl . Therefore, for a family of purities (jTEI), it is possible to conclude, that 
"superpurity" and entropy-based purity lead to limiting cases of all multidimensional 
uncertainty relations. 

Appendix C. 

Starting from the equation ^K)^ . with a sequence of nonnegative numbers and M an 
arbitrary real nonnegative constant, we write 



m=0 



2M + n 



2 J2 U{M -m) + Y^ em(2m + 



n] 



m=0 



771=0 



The following (in)equalities hold: 



m] 



1. (C.l) 



(C.2) 



m=0 



0<m<M 



J2 gt\M-m)U9. 



in) 
m 



0<m<M 



0<m<M 



J2 9t\M-m 

00 



1/r 



1 Vp 



.O^m<A/ 



(C.3) 



(C.4) 



m=0 
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with two real parameters p,r ^ 1, 1/p + 1/r = 1. The equahty sign in relations 
()U.2|) and ()C.4|1 holds if = for m > M. Relation ()C.3|) changes to equality, when 
Cm oc gm\M — m)''"^ in the interval ^ m ^ M. (In)equality (jdSj) is a general form 
of the Holder inequality for the weighted sum [30], see also [31]. 
Combining the (in)equalities (jCy.ljl - ()(14j) gives a relation 



2M + n 



2B{M, n, r)^/>g + ^ Ui'^m + 



n) 



where 



B{M,n,r)= J] g^^\M 



m=0 



m] 



> 1, 



0^m<M 



m=0 



1/rp 



(C.5) 



(C.6) 



(C.7) 



From the condition 1/p + 1/r = 1 it follows that p = r/(r — 1) and then the 
[in)equality (PHj) results. 



Appendix D. 

In order to find an integral in approximation (j30|) . we start from introduction of a new 
variable x = M — m 

1 

B(M,n,r)^- -/ dx (M - (D.l) 

[n - 1)! Jo 

Application of binomial formula to (M— x)"~^ and interchanging the order of integration 
and summation leads to 

B(M.„.r)«M— 'g^^— iZ^i^-^ (D.2) 

The last sum can be calculated by use of formula (5.41) from Graham, Knuth and 
Patashnik book [32], leading at last to ()3H). 
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